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We take a two-step theoretical approach to study magnetism of rare earth quasicrystals by con¬ 
sidering Ising spins on quasiperiodic tilings, coupled via RKKY interactions. First, we compute 
RKKY interactions from a tight-binding Hamiltonian defined on the two-dimensional quasiperiodic 
tilings. We find that the magnetic interactions are frustrated and strongly dependent on the local 
environment. This results in the formation of clusters with strong bonds at certain patterns of the 
tilings that repeat quasiperiodically. Second, we examine the statistical mechanics of Ising spins 
with these RKKY interactions, using extensive Monte Carlo simulations. Although models that 
have frustrated interactions and lack translational invariance might be expected to display spin 
glass behaviour, we show that the spin system has a phase transition to low-temperature states 
with long-range quasiperiodic magnetic order. Additionally, we find that in some of the systems 
spin clusters can fluctuate much below the ordering temperature. 

PACS numbers: 75.50.Kj, 71.15.-m 


Quasicrystals with their unusual atomic structure 
characterized by long-range order without a three- 
dimensional translational periodicity^ are known to have 
rather exotic physical properties, and so far we lack a 
good theoretical understanding for many of them. For 
instance, the electronic properties of this material class 
include a pseudogap at the Fermi energy^i^ and multi¬ 
fractal (critical) wave functions^ resulting in anomalous 
electronic transport^. An important question in this 
context is how these electronic properties influence the 
magnetism in quasicrystals 

In general, one can distinguish two types of magnetic 
quasicrystals: those with magnetic moments at (i) tran¬ 
sition metal or (ii) rare earth sites. In the first, moment 
formation is an important aspect of the theoretical prob¬ 
lem since they appear only at a small fraction of sites.— By 
contrast, the second class suggests a simpler description, 
with well-defined local moments at concentrations of 5- 
10% interacting via long-range RKKY interactions which 
are mediated by the conduction electrons— Examples 
of the latter class include the icosahedral i-ZnMgR and 
i-AgInR materials^iiii, as well as decagonal d-ZnMgR 
materialsi^ and the recently discovered binary phase&i^ 
i-RCd. 

Here we use tight-binding models and Ising models 
to address this theoretical problem. We examine the 
form of the RKKY interactions based on a tight-binding 
Hamiltonian defined on two-dimensional quasiperiodic 
tilings (see Section HI. We find that the coupling between 
pairs of sites depends not only on their distance but also 
varies strongly with the local environment on the tiling. 
Although we find ferromagnetic and antiferromagnetic 
bonds as in periodic systems, the magnetic interactions 
do not show a well-defined spatial period with a Fermi 
wave vector because quasicrystals do not have a Brillouin 
zone. Second, we study the consequences of these inter¬ 
actions (see Section |H]), by taking them as exchange con¬ 
stants between Ising spins located at a fraction of sites 
on the tiling (those with a particular coordination num¬ 


ber). As all systems show a combination of frustration 
and aperiodicity one might expect to observe a spin glass 
at low temperatures!^, and in other settings quasiperi¬ 
odic systems are known to behave like random onesi^. 
Our results exclude canonical spin glass behaviour via 
the temperature dependence of the order parameter sus¬ 
ceptibility, the overlap susceptibility, and the magnetic 
structure factor of the ground state. Instead we find 
that the ground state consists of repeating small clusters 
of spins, with strong interactions within each cluster and 
weaker couplings between different clusters. Analysing 
the magnetic structure factor we find that these systems 
show a transition to a state with long-range magnetic 
order at low temperatures. 

While there has been extensive previous work study¬ 
ing spin models for magnetic quasicrystals, to the best 
of our knowledge none has used RKKY interactions 
computed from a quasiperiodic electron Hamiltonian. 
Moreover much of this previous work omits the frus¬ 
tration effects that are a natural consequence of long- 
range oscillatory RKKY interactions. Nearest neigh¬ 
bour exchange on a bipartite tiling necessarily leads 
in an Ising antiferromagnet— ^ to a classical ground 
state with two-sublattice order. In Heisenberg mod¬ 
els it is known that this order may survive quantum 
fluctuations!^ or the inclusion of dipolar and further 
neighbour interactions— RKKY interactions with 
a form taken from periodic systems have been used in 
studies of Ising models for the Penrose tiling^S and for 
i-ZnMgHo models^i2^, yielding antiferromagnetic or fer- 
rimagnetic ordered ground states. Neither of these mod¬ 
els incorporates the expected unique coupling of the ex¬ 
otic electronic properties to the magnetism in quasicrys¬ 
tals. Our model is designed to address this aspect of the 
physics in a simple way. 

In previous work (Ref. [H) we have studied the struc¬ 
ture of the RKKY interactions for a limited number of 
systems. Here we present a much more thorough analy¬ 
sis and classification of quasiperiodic spin systems. An 
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important improvement is a new algorithm for the com¬ 
putation of the RKKY interactions which has higher ac¬ 
curacy and is also applicable to systems with gaps in 
the density of states, in contrast to the previously em¬ 
ployed continued fraction expansion. We use the new 
algorithm to compute the RKKY interactions for a large 
set of systems. This allows us to identify different classes 
of behaviour with respect to the structure of the mag¬ 
netic ground state. The new results also provide a more 
detailed insight into the nature of low-temperature fluc¬ 
tuations. Finally, we apply finite size scaling to show 
that the critical behaviour is consistent with the two- 
dimensional Ising universality class. 


I. RKKY INTERACTIONS IN TIGHT-BINDING 
MODELS 

Introducing Ising spins crj at some of the sites I of 
the tiling, the spin Hamiltonian has the form FIrkky = 
Here, A represents the coupling of the local 
moment to the conduction electrons. We model the con¬ 
duction electrons by the tight-binding Hamiltonian Hei = 
m) 10 (■'^1 with one orbital per site and equal hopping 
amplitudes between all nearest neighbours of a quasiperi- 
odic tiling (see Fig. [T|). The local susceptibility^i for 
T = 0 and Fermi energy E-p is (for a review see Ref. [25ll 

Xi,m = - f O' [Gi^^G^^i] sign(£' - £'F)d£;. (1) 

J — OO 

Here Gi^rn = (^| G \m) is a matrix element of the retarded 
Green function G = [E + iQ — for the conduction 

electrons. Writing the Green functions in terms of the 
eigenstates as 


^ E-Ea+lQ 

a. 


( 2 ) 


Eq. ([T]) evaluates to 


Xl,m 


2y 


Ea) , 


( 3 ) 

where the sum runs over all eigenstates with either E^ < 
Ep < Ep or Ep < Ep < Ea- This equation clearly 
shows that the properties of the electronic eigenstates of 
the Hamiltonian i?ei have an important influence on the 
RKKY interactions. In the case of quasicrystals, this ap¬ 
proach also naturally includes the multifractal structure 
of the electronic wave functions 

The numerical evaluation of Eq. m leads to a di¬ 
vergence of the susceptibility if the Fermi energy Ep 
is close to an energy eigenvalue Ea- To overcome this 
problem we approximate the sign-function in Eq. [T] with 
tanh((i? — Ep)/ri), where rj can be interpreted as a tem¬ 
perature parameter. This yields a good approximation 
as long as the band width of the system ^ rj. We use 
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FIG. 1. Tilings and RKKY interactions Xi.m for magnetic mo¬ 
ments at sites with coordination number z, coupling strength 
A = 1, and rj = 0.03: (a) 5th approximant of the Pen¬ 

rose tiling, showing only interactions above the threshold 
\Xi,m\ > 0.005, z = 5, and Ep = 0.07; (b) 4th approximant of 
the Ammann Beenker tiling with |xi.m| > 0.006, 2 = 4, and 
Ep = 1.12. 


r] = 0.03 for all calculations which is much smaller than 
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the band width for the considered quasiperiodic systems. 
Re-evaluating Eq. [TJ we can then rewrite Eq. [3] replacing 
the sign(E'F — Ea) with the expression 

sinh ( 

_ V ^ ^ _ (4) 

2 cosh cosh 

Our approach is based on the direct diagonalization 
of the electron Hamiltonian He\. Using fast numerical 
libraries (lapack) and parallelization we are able to com¬ 
pute the magnetic interactions for systems with up to 10^ 
sites. Our method is more accurate than the continued 
fraction expansion, which has been employed previously 
for quasiperiodic tilings 

The method can be applied to general tight-binding 
systems. For a particular model we need to specify three 
features: (i) the tiling; (ii) the positions of the spins on 
the tiling; and (iii) the value of the Fermi energy Ep. In 
this paper we study two different quasiperiodic systems. 

(A) Penrose tiling (see Fig.[T^): This tilings can be used 
to describe the structure of decagonal quasicrystals.— It 
has 5 different local environments with 3 to 7 nearest 
neighbours. We consider the 4th to 6th approximants 
with 521, 1364, and 3571 sites respectively. 

(B) Ammann Beenker tiling (see Fig. [Ud): This tiling 
models the structure of octagonal quasicrystals^ and 
possesses 6 different local environments with 3 to 8 near¬ 
est neighbours. We consider the 4th and 5th approximant 
with 1393 and 8119 sites respectively. 

As stable quasicrystals are characterized by a high degree 
of order—, we choose the spins to be at sites with a fixed 
coordination number z. 

An example for the susceptibility Xi,'m is shown in 
Fig. m for each of the two tilings. In general, the suscep¬ 
tibility strongly depends on Ep and on the local environ¬ 
ment, and interactions of similar strength can be found 
whenever the same local pattern reoccurs in the tiling. 
Since any local pattern of linear dimension L is repeated 
in a distance 0{L) for each quasiperiodic tiling)^ the 
strongest bonds form clusters that have a fixed form and 
are repeated quasiperiodically. 

In our previous work (Ref. 12311 we studied the distribu¬ 
tion of Xi,m in detail for the Ammann Beenker tiling. We 
observe the same qualitative behaviour also for the Pen¬ 
rose tiling and we briefly summarize the main findings 
here before focussing on results from Monte Carlo simu¬ 
lations. In general, the susceptibility xi.m is proportional 
to the local DOS at Ep which is known to vary strongly 
with the local environment^. This also means that inter¬ 
actions are smaller when Ep lies in a pseudogap which is a 
typical feature of quasicrystals. We find oscillations be¬ 
tween ferromagnetic and antiferromagnetic interactions 
as a function of site separation r, measured along the 
shortest bond path. This behaviour differs significantly 
from that in periodic systems with a spherical Fermi sur¬ 
face, where interactions oscillate within a power-law en¬ 
velope that varies as r~'^. Although in some cases the 
average (xi.m) is reasonably well described by this power 


law, individual interactions in quasiperiodic tilings can 
be considerably larger, and we usually observe a wide 
range of interaction strengths at each r. This results in a 
quasi-random contribution to the magnetic interactions 
due to the many local environments in a quasicrystal. 


II. MONTE CARLO SIMULATION 


We use Monte Carlo simulations to study the statisti¬ 
cal mechanics of the spin Hamiltonian TIrkky • We apply 
parallel tempering which simulates multiple copies of the 
system at different temperatures T^, so reducing corre¬ 
lation times— We use a geometric distribution of the 
temperatures at low temperatures with Tm = c™ Tmin 
and a uniform spacing with Tm = Tm-i + C 2 above the 
phase transition. For each replica one Monte Carlo Sweep 
consists of N single-spin flip attempts and one replica- 
swap attempt, where N is the number of spins in the 
system. The number of replicas (typically 60-100) is cho¬ 
sen to have an acceptance ratio of at least 30% for the 
replica swaps. 

We compute a selection of observables to study the 
phase transition and the ground state of this system. 
This includes the energy E = ^Xi,mO'i'^m and the 

magnetization M — Thermal averages are com¬ 

puted over 5 x 10^ Monte Carlo Sweeps. The heat capac¬ 
ity per spin is C = [{E"^) — (E)"^] and the suscepti¬ 
bility per spin is y = [(Af^) — (|A/|)^]. The former 

provides a signal for phase transitions. While for conven¬ 
tional antiferromagnets the staggered magnetization is a 
useful tool to detect long-range order, we do not know the 
structure of the ground state a priori. To look for long- 
range magnetic order in an unbiased way, we first search 
for the ground state configuration {^i} and then compute 
the corresponding order parameter Mgs = 
can use M™ to define the order parameter susceptibility 
Xop = ^I(Mp-(|Mgs|)2]. 

Further, we consider some quantities commonly used 
to study spin glasses which turn out to be sensitive to 
different ordering patterns. This includes the overlap 
q = ^ it® distribution P{q) obtained from 

the simulations of two independent replicas with the same 
interactions and temperature. We find that the Binder 
cumulant Bsg = ^(3 ~ (<Z‘^)/('Z^)^) is a good tool to iden¬ 
tify order, whether it is (anti)ferromagnetic or spin-glass¬ 
like. This becomes clear by looking at the limiting cases: 
at high temperatures the spins are uncorrelated and we 
obtain q ^ and RsG = 0; in contrast, if there 

is only a pair of ground states related by a global spin 
inversion at low temperatures, we obtain q = ±1 and 
Bsg = 1 - 
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TABLE I. Classification of low-temperature states for differ¬ 
ent Fermi energies Ep and spins at sites with coordination 
numbers 2 for the Ammann Beenker tiling and the Penrose 
tiling. (We did not study the Penrose tiling with 2 = 4 be¬ 
cause this system has a very inhomogeneous distribution of 
spins.) Abbreviations: FMR - ferromagnetic regions, AFl - 
antiferromagnetic with a single ground state, AF2 - antifer¬ 
romagnetic with multiple ground states. 


system 

Ef 

z = 3 

2=4 

2 = 5 

Penrose 

-4 

FMR 

- 

FMR 


-3.5 

AFl 

- 

AF2 


0 

AFl 

- 

AFl 


0.07 

AF2 

- 

AFl 


2.33 

AF2 

- 

AFl 


2.79 

AFl 

- 

AF2 

Ammann 

-4 

FMR 

FMR 

AFl 

Beenker 

-3.5 

AF2 

AF2 

AFl 


0 

AFl 

FMR 

AFl 


0.7 

AF2 

AF2 

AFl 


1.12 

AFl 

AF2 

AF2 


1.3 

AF2 

AF2 

AF2 


1.95 

AF2 

AFl 

AF2 


2.7 

AF2 

AFl 

AF2 


A. Phase Transition 


We use Monte Carlo simulations for different parame¬ 
ter sets for the Ammann Beenker tiling and the Penrose 
tiling obtained by varying Ep and choices 2 for the co¬ 
ordination number of the magnetic sites. Our results 
are briefly summarized in Table HI We were able to clas¬ 
sify the systems by the structure of their ground state 
into two main classes: The first class is characterized 
by ground states with large ferromagnetic regions and 
we denote it as class FMR. The second class (class AF) 
shows ground states with antiferromagnetic correlation. 
We can further subdivide this class with respect to the 
number of ground state. While class AFl has only a sin¬ 
gle pair of ground states related by time-reversal symme¬ 
try, the class AF2 is characterized by multiple low-energy 
states. This leads to distinct differences in the statisti¬ 
cal mechanics of the systems as outlined below. In the 
following we discuss the classes in detail. 

Due to the computational complexity of our simula¬ 
tions, we can only obtain simulations results for very 
few successive approximants for each system (see also 
Sec. Ill Dl) . For these successive approximants we find the 
same qualitative behaviour provided they are big enough 
(about 300 spins) to capture the physics of the quasiperi- 
odic system. 


1. FMR class 

Systems with large ferromagnetic regions mainly oc¬ 
curs for Fermi energies near the band edges. In these 
cases, the dominant contributions to xi,m are from elec¬ 
tronic wave functions which have a slowly varying enve¬ 
lope and components on the two sublattices of the tiling 
with either the same sign (upper band edge) or opposite 
signs (lower band edge). Both cases lead to dominant 
ferromagnetic interactions between the spins with very 
weak frustration. 

A typical example is the Penrose tiling illustrated in 
Fig- [2] for the Fermi energy Ep = —4 with N = 693 
Ising spins at sites with coordination number z = 3. The 
ground state shows large ferromagnetic regions and the 
weak long-range antiferromagnetic bonds lead in this case 
to the formation of large oppositely-aligned regions. Of¬ 
ten the spin patterns form rather elongated regions in 
order to avoid domain walls along strong ferromagnetic 
bonds. According to the results in Fig. [5^, the system 
shows a phase transition at about Tf Ri 0.02 from a para¬ 
magnetic state at high temperatures to an ordered state 
at low temperatures, as suggested by the broad peaks in 
the heat capacity C and the susceptibility y as well as 
an increase of the Binder cumulant. At very low tem¬ 
peratures the susceptibility y and the order parameter 
susceptibility yop show a 1/T-divergence due to the fluc¬ 
tuations of small clusters or single spins. This is discussed 
in more detail in Sec. EC] 


£. AF classes 

Most of the studied systems belong to the AF class 
for which we find a phase transition to low-temperature 
states with antiferromagnetic correlations. In the follow¬ 
ing we focus on this class, and in Sec. Ill Bl we show that 
the low-temperature state has quasiperiodic Neel order. 
To illustrate this transition, we first present simulation 
results for both sub-classes, AFl and AF2. 

A typical example for the AFl class is shown in Fig. [31 
It corresponds to the Penrose tiling in Fig. [TJi with N = 
375 Ising spins at sites with coordination number 2 = 5 
and taking Ep — 0.07. For this system, the heat capacity 
C has a peak near Tf r: 0.05, and the Binder cumulant 
Bsg approaches 1 below this temperature. This suggests 
that the system is paramagnetic for T > Tf, and that 
at Tf a macroscopic fraction of spins lock together into a 
state with long-range rigidity. There is a large peak in the 
order parameter susceptibility yop at Tf associated with 
the spontaneous breaking of the global Ising symmetry of 
the spin model when spins order into a low-temperature 
phase with nearly zero magnetization M but non-zero 
Mgs. In our previous work (see Ref. [23l we found the 
same qualitative features for the Ammann Beenker tiling 
with Ep = 1.95 and 2 = 4. 

In contrast, behaviours in the AF2 class are more com¬ 
plex. An example is the Ammann Beenker tiling with 
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temperature r temperature T 



-0.004 -0.002 0.000 

susceptibility x,.„ 



FIG. 2. Simulation results for the Penrose tiling with 2 = 
3 and Ey = —4 (class FMR): (a) different observables and 
(b) lowest energy spin configuration (red and blue circles for 
up and down spins) and (c) overlap distribution Pdql) at 
different temperatures. 


N = 478 spins at sites with z = 4 and Ey = 1.12 in 
Fig. [T]d. The temperature behaviour of the observables 
and ground state are shown in Fig. SI Clear indications 
for a phase transitions come from the Binder cumulant 
SsG which increases from 0 to values above 0.95 below 
Tf Ri 0.04. Again this is due to the transition from a 
paramagnetic state for T > Tf to a state with broken 



FIG. 3. Simulation results for the Penrose tiling with z = 5 
and Tf = 0.07 (class AFl): (a) different observables and (b) 
lowest energy spin configuration (red and blue circles for up 
and down spins) with strongest bonds according to Fig. [1^ 
and (c) overlap distribution Pd?!) at different temperatures. 


Ising symmetry with antiferromagnetic correlations be¬ 
low Tf. The heat capacity C and the susceptibility also 
show very broad maxima near Tf. 

To distinguish between the two AF subclasses, we show 
the overlap distribution P{q) in Figs. [51: and 01;. For 
both systems P{q) is centred at g = 0 for T > Tf but 
develops one or more pairs of peaks at ±qi{T) for T < 
Tf. If we observe a single pair of peaks go —>■ ±1 at 
very low temperatures, the spins freeze in a particular 
ground state. The system in Fig. [ 51 : is an example of 
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temperature T 



FIG. 4. Simulation results for the Ammann Beenker tiling 
with 2 = 4 and Ey = 1.12 (class AF2): (a) different observ¬ 
ables and (b) lowest energy spin configuration (red and blue 
circles for up and down spins) with strongest bonds accord¬ 
ing to Fig. [Ha and (c) overlap distribution P{\q\) at different 
temperatures. 


this behaviour. In contrast, the overlap distribution in 
Fig. |4]: shows 4 peaks at qq = ±1 and q ~ ±0.75 for 
T —^ 0, indicating that there are two competing pairs of 
low-temperature states. We use the standard deviation 
of the overlap distribution a{P{\q\)) at low temperatures 
to distinguish the two classes. We find that a — 0.02 is 
a good threshold value to distinguish the AF2 class with 
multiple low-temperature states (cr > 0.02) from the AFl 
class with a single ground state (cr < 0.02). 


For the AF2 class, signatures of the multiple ground 
states are also clearly visible in the order parameter sus¬ 
ceptibility Xop which diverges at low-temperature due to 
a Curie-like contribution originating from the fluctua¬ 
tions of the system between the different ground states 
(see Fig. SI). Also the susceptibility x of this system 
has an additional rather sharp peak at Ti r; 0.005. In 
Sec. Ill Cl we show that this peak is caused by fluctuations 
of loosely-coupled spins. 

Note that the overlap distributions we obtain are strik¬ 
ingly different from those observed in spin glass simu¬ 
lations using the three-dimensional Edwards Anderson 
model or the Sherrington-Kirkpatrick model^ for which 
P(q) has much more weight at smaller q below Tf. 

Although we can access only a few approximants nu¬ 
merically, we expect that the number of low-temperature 
states in AF stays finite in the thermodynamic limit. Ev¬ 
idence for this comes from the finite size scaling results in 
Sec. lIIDi that show an excellent agreement of the different 
approximants according to relevant scaling equations. 

B. Structure of Ground State 

A key question is what type of order characterises the 
low-temperature phase. Short-distance correlations are 
readily apparent in the ground state spin configuration. 
Eor example, the Penrose tiling in Fig. [3}) contains many 
pentagons with a ferromagnetic spin configuration form¬ 
ing on the 5-fold symmetric stars of the tiling. Long- 
range order results from a tendency of neighbouring clus¬ 
ters to anti-align due to antiferromagnetic bonds between 
the clusters. For the Ammann Beenker tiling in Fig. |41 d 
two distinct local patterns can be identified: octagons 
with an antiferromagnetic spin configuration form on the 
8-fold star patterns, and small V-shaped clusters form 
along the strong ferromagnetic bonds. Again long-range 
order is caused by a tendency of neighbouring clusters to 
anti-align. 

In both cases the nature of this order is revealed 
to be quasiperiodic by the magnetic structure factor 

computed from the spin con¬ 
figuration of the ground state. In contrast to pe¬ 
riodic systems, quasicrystals do not posses a Brillouin 
zone and understanding the diffraction pattern is usu¬ 
ally challenging. However, it is possible to derive the 
nuclear structure factor of quasiperiodic tilings theoret¬ 
ically, and one can show that the quasiperiodicity leads 
to Bragg peaks in a pattern that forms a reciprocal- 
space quasiperiodic tilingi^ The corresponding length 
scale is 47rT/(56) for the reciprocal-space Penrose tiling^S 
and 1 / (2b) for the reciprocal-space Ammann Beenker 
tiling^i^ in units of the real-space bond length b and 
the golden mean t = (1 ± •\/5)/2. When only a subset of 
sites of the tiling is chosen, the nuclear structure factor 
is rather robust and there are often only minor changes 
in the amplitudes of the Bragg peaks 

Numerical results for a ferromagnetic state (all = 1) 
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on the Penrose tiling with z = 5 are shown in Fig. [5^. 
This gives the nuclear structure factor, and the diffrac¬ 
tion intensity M(k) indeed shows a pattern of peaks 
which can be overlaid by a Penrose tiling with the ex¬ 
pected length scale. All sites of the reciprocal-space Pen¬ 
rose tiling fall on one of two different patterns: a high- 
intensity peak or the centre of a ring with 10 medium- 
intensity peaks as highlighted in Fig. [S^. The 10-fold 
rotational symmetry of the Penrose quasicrystal is also 
visible around the central Bragg peak at k = 0. 

In the case of an antiferromagnetic ground state the 
correlations gives rise to selection rules that result in a 
shift of the reciprocal space indices.— For the example 
in Fig. [5}), this means that all peaks split into several 
peaks along the principal directions of the tiling. This 
is clearly visible for the central Bragg peak which splits 
into 10 sizeable peaks along the principle directions of 
the Penrose tiling. All other high-intensity peaks for the 
ferromagnetic state in Fig. [5^ also split into 10-fold rings 
of medium-intensity peaks for the ground state with an¬ 
tiferromagnetic order. The opposite effect occurs for the 
10-fold rings with medium-intensity peaks of the ferro¬ 
magnetic state. While each of these peaks again splits 
into 10 peaks along the principal directions, constructive 
interference leads to the formation of a high-intensity 
peak in the centre of the ring in the antiferromagnetic 
ground state. The annotations in Fig. [5] highlight this 
transformation of the diffraction patterns (i.e. a peak to 
a 10-fold ring and vice versa) between the ferromagnetic 
and antiferromagnetic spin configuration. The observed 
structure factors are also consistent for different approx- 
imants, i.e., the peaks occur at the same positions and 
become sharper for bigger approximants. 

A possible concern is whether the quasiperiodic site 
locations by themselves might be sufficient to generate 
these features in the diffraction pattern, without mag¬ 
netic order. To test this we have computed the magnetic 
structure factor for a random spin configuration on the 
tiling (Fig. [St): its amplitude fluctuations are Gaussian^^ 
and it does not show any of the high intensity peaks that 
are present for the ground state and for a ferromagnetic 
one. This clearly indicates that the observed pattern is 
due to the long-range magnetic order and not due to the 
atomic order. 


C. Formation and Fluctuations of Spin Clusters 

Our results show some distinctive features in addition 
to those arising from ordering. In all systems considered, 
the RKKY interactions generate strongly coupled spin 
clusters with weaker inter-cluster couplings as illustrated 
in Fig. dl For a given choice of magnetic site and Ap 
the most prominent clusters have a fixed form that is 
frequently repeated. According to Conway’s theorem, 
the prominent clusters repeat quasiperiodically because 
any local pattern of linear dimension L is repeated in a 
distance 0(L) for each quasiperiodic tiling— 





FIG. 5. Magnetic structure factor (scaled as M(k)i for better 
visibility) of the 6th approximant of the Penrose tiling with 
N = 1003 spins on sites with coordination number z = 5: 
(a) ferromagnetic state, (b) lowest energy spin configuration 
for Ep = 0.07, and (c) random spin configuration. In the 
first two cases a Penrose tiling with the expected length scale 
47rr/(56) ~ 4.07 for 6 = 1 is overlaid on the right half of the 
image to visualize the quasiperiodic structure. Note that the 
peaks and 10-fold rings (marked by arrows) are exchanged 
between (a) and (b). 
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FIG. 6. Simulation results for the Ammann Beenker tiling 
with 2 = 3 and Ey = 1.12 (class AFl): (a) different ob¬ 
servables, (b) lowest energy spin configuration (red and blue 
circles for up and down spins), and (c) overlap distribution. 


In addition, it is striking that in many systems some 
spins or small clusters are free to fluctuate at temper¬ 
atures much below the ordering transition, yielding a 
Curie-like contribution to % and/or Xop- A typical ex¬ 
ample is the Ammann Beenker tiling with z = 3 and 
E-p = 1.12 shown in Fig. [51 The system belongs to 
the class AFl with a single ground state. At about 
Tf Ri 0.06 the system shows a (quite broad) phase transi¬ 


tion to a state with antiferromagnetic order. In contrast 
to the previous systems, some spins are able to fluctuate 
much below the ordering temperature resulting in a \/T- 
divergence of the susceptibility %. We expect that weak 
interactions will suppress spin fluctuations at low tem¬ 
peratures, so that X approaches zero as T —>■ 0. Hence, 
in some systems x shows a peak at low temperatures if 
spin fluctuations are present. 

To study further these fluctuations we compute for 
each spin i the thermal average of the spin-direction 
with respect to the ground state using (rUgs) = 
(sign(Mgs) For a paramagnetic state at T > Tf 

with no preferred spin direction we expect (rngg) = 0. In 
contrast, we get (rUgg) = 1 when the system freezes into 
a single ground state. We find that many of the studied 
systems have some ’loose spins’ that fluctuate strongly 
even for T <C Tf. In Fig. [7] we show the spin fluctuations 
(Wgg) for the Ammann Beenker tiling in Fig.[S]at different 
temperatures. We find that even at very low tempera¬ 
tures single spins on many of the 8-fold rings and some of 
the semi-circular clusters are able to fluctuate. Looking 
at the magnetic interactions this is caused by frustration 
effects due to presence of ferromagnetic and antiferro¬ 
magnetic interactions and the rather broad distribution 
of interaction strengths with only weak inter-cluster in¬ 
teractions. 

Fluctuations due to loose spins are also present in the 
Ammann Beenker tiling shown in Fig. |T1 as can be seen 
from the peak in the susceptibility x at low-temperatures. 
However, this system shows a second type of fluctuations, 
between the two almost degenerate ground states. This 
leads to the divergence of the order parameter suscepti¬ 
bility Xop at low temperatures as can be seen in Fig. 2^. 
This is an artefact of the parallel tempering method 
which allows fluctuations between different ground states 
that are not possible with realistic systems. This effect 
can be observed for all systems in the AF2 class. 


D. Finite Size Scaling 

Finite size scaling is a valuable tool for studying phase 
transitions in finite systems and extrapolating to the 
thermodynamic limit. However, it is difficult to employ 
it for quasicrystals as on the one hand the system size of 
successive periodic approximants grows fast and on the 
other hand we require fairly large system sizes to see the 
effects of the quasiperiodicity. For instance, the size of 
successive approximants grows with a factor t = 3 + 2\/2 
for the Ammann Beenker tiling and r = (1 -|- •\/5)/2 for 
the Penrose tiling. Further complication are caused by 
the fluctuations of loose spin, multiple ground states and 
rather broad transitions in some systems, as all of them 
lead to a complex temperature behaviour of the relevant 
observables. 

Luckily, some systems avoid these problems. One ex¬ 
ample is the Penrose tiling in Fig. 2] with Ep = 0.07 
and z = 5 for which we were able to compute the ob- 
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FIG. 7. Fluctuations (fUgg) of spins for the system in Fig. |6] at different temperatures: (left) for T = 0.01 loose spins and 
semi-circles of spins are able to fluctuate leading to the peak in x\ (middle) at T = 0.031 additional fluctuations occur, and 
(right) at the transition temperature T = 0.078 < T{ fluctuations are present at all sites but {rn\s) > 0. 




scaled temperature L{T-T^ 

FIG. 8. Finite size scaling results of the order parameter sus¬ 
ceptibility Xop and the Binder cumulant Bsg for the Penrose 
tiling with Ef = 0.07 and 2 = 5. Results are shown for the 
5th approximant with N = 375 and the 6th approximant with 
N = 1003 spins. 


servables for the 5th and 6 th approximants, with 375 
and 1003 spins respectively. We also obtained results 
for smaller approximants but they show quite different 
transition temperatures because they are too small to 
adequately describe the quasiperiodic structure. For the 
two large approximants, we find that the peak of the 
order parameter susceptibility Xop becomes sharper and 
the slope of the Binder cumulant B^g increases with the 
system size. 

While in an infinite volume the correlation length ^ 
diverges near the transition point as ^ oc with 


t = T — Ti, in simulations with finite size N = L‘^ the 
system already becomes effectively ordered for correla¬ 
tion lengths ^ « L. Therefore, the observables can be 
described by the finite size scaling equations XopiL, T) = 
for the order parameter susceptibility 
and Bsg{L,T) = gB{L^^''\t\) for the spin glass Binder 
cumulant We expect our systems to be in the 2D 
Ising universality class for which the exact exponents are 
known to be jz = 1 and 7 = 7/44^ To determine an es¬ 
timate for the critical temperature Tf we use the Binder 
cumulant Bsg which is less sensitive to finite size effects. 
For the Penrose tiling with Ep = 0.07 and 2 = 5, we find 
an excellent agreement of the scaled Binder cumulant for 
the two approximants with Tf = 0.0445 (see Fig.|S]). Scal¬ 
ing of the order parameter susceptibility is also in very 
good agreement for this critical temperature. 


E. Comparison to Experimental Results 

So far quasiperiodic magnetic order has not been ob¬ 
served in rare earth quasicrystals 4 Experimental data 
for many of these materials shows a spin-glass-behaviour 
with a separation of the field-cooled and zero-field-cooled 
magnetic susceptibility below a freezing temperature 
Tf44J4 In most cases the Weiss temperature 0 is neg¬ 
ative, indicating predominantly antiferromagnetic inter¬ 
actions. Also —Q/Tf is rather large with values of 5 to 
10 implying strong frustration4d4 

However, some 1/1 cubic approximants are known to 
have long-range magnetic order. For instance, the CdgR 
approximants show long-range antiferromagnetic order 4^ 
Interestingly, the 1/1 cubic approximant of AuSiGd and 
AuGeGd are both ferromagnetic below a freezing tem¬ 
perature but at even lower temperatures the AuGeGd 
approximant alone shows an additional transition to a 
spin-glass4^ According to current structure models, the 
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major difference is some site disorder of the rare earth 
atoms in AuGeGd which is not present in the AuSiGd 
approximant.— In general, site disorder of the rare earth 
atoms is very common in quasicrystals. For instance, 
the structure model for the recently discovered binary 
quasicrystals GdR also includes site disorder for some of 
the rare earth positions . This suggests that the spin 
freezing observed in rare earth quasicrystals may be at¬ 
tributed to structural disorder. 

With respect to observed spin freezing in rare earth 
quasicrystals, some features are in agreement with canon¬ 
ical spin glasses but there are also important differences. 
The frequency dependency of the a.c. magnetic suscep¬ 
tibility and the behaviour of the non-linear suscepti¬ 
bility in i-ZnMgRe are typical for a conventional spin- 
glass transition.® However, neutron scattering data for 
i-ZnMgTb and i-ZnMgHo show the formation of short- 
range order within spin-clusters at about 20K. Although 
individual spins are frozen within the clusters, the clus¬ 
ters are able to fluctuate in a rather broad temperature 
range of about 5.8K < T < 2GK^— The latter ob¬ 
servation fits well with our results for the formation of 
strongly coupled spin clusters on certain patterns of the 
tiling and the fluctuations of these clusters much below 
the ordering transition being permitted by the weaker 
inter-cluster interaction. 

A natural extension of our model is to include site dis¬ 
order in the systems and study their magnetic proper¬ 
ties. We find in preliminary work that RKKY interac¬ 
tions computed in our model show a high sensitivity to 
a weak random on-site potential. This also implies that 
even disorder restricted to non-magnetic sites will influ¬ 
ence the couplings between spins. 

III. CONCLUSION 

We have used Ising spins with RKKY interactions com¬ 
puted from tight-binding models on quasiperiodic tilings 


as a simple caricature for the magnetic properties of rare- 
earth quasicrystals. The RKKY interactions were com¬ 
puted with an improved numerical method which directly 
uses the multifractal electronic eigenstates of the system. 
For all systems we find the emergence of strongly-coupled 
spin clusters with weak inter-cluster coupling on certain 
patterns of the tiling which repeat quasiperiodically. 

Using extensive Monte Carlo simulations, we find that 
all studied systems show a phase transition at a temper¬ 
ature Tf at which the global Ising symmetry of the spin 
model is broken spontaneously. Evidence for this tran¬ 
sition includes the temperature-dependence of the order 
parameter susceptibility, spin glass Binder cumulant, and 
the overlap distribution. While some systems have large 
ferromagnetic regions at low temperatures, the major¬ 
ity of systems show a transition to a quasiperiodic Neel 
state. We demonstrate that the low-temperature state 
has long-range quasiperiodic magnetic order by analysing 
the magnetic structure factor. This result is striking tak¬ 
ing into account the frustration and quasi-randomness of 
the magnetic interactions in our quasicrystal model. 

The formation of strongly coupled clusters and their 
fluctuations even at very low temperatures appears to be 
consistent with experimental observations of the forma¬ 
tion of ordered spin clusters and collective spin fluctua¬ 
tions. In contrast, the nature of the ordering transition 
in the model shows clear differences to experiments, and 
it remains an interesting open problem to investigate the 
origin of the spin freezing found in many of these mate¬ 
rials. 
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